Population-specific Mutation Patterns in Breast Tumors from African American, European American, and Kenyan Patients.

Women of African descent have the highest breast cancer mortality in the United States and are more likely than women from other population groups to develop an aggressive disease. It remains uncertain to what extent breast cancer in Africa is reminiscent of breast cancer in African American or European American patients. Here, we performed whole-exome sequencing of genomic DNA from 191 breast tumor and non-cancerous adjacent tissue pairs obtained from 97 African American, 69 European American, 2 Asian American, and 23 Kenyan patients. Our analysis of the sequencing data revealed an elevated tumor mutational burden in both Kenyan and African American patients, when compared with European American patients. TP53 mutations were most prevalent, particularly in African American patients, followed by PIK3CA mutations, which showed similar frequencies in European American, African American, and the Kenyan patients. Mutations targeting TBX3 were confined to European Americans and those targeting the FBXW7 tumor suppressor to African American patients whereas mutations in the ARID1A gene that are known to confer resistance to endocrine therapy were distinctively enriched among Kenyan patients. A Kyoto Encyclopedia of Genes and Genomes pathway analysis could link FBXW7 mutations to an increased mitochondrial oxidative phosphorylation capacity in tumors carrying these mutations. Finally, Catalogue of Somatic Mutations in Cancer (COSMIC) mutational signatures in tumors correlated with the occurrence of driver mutations, immune cell profiles, and neighborhood deprivation with associations ranging from being mostly modest to occasionally robust. To conclude, we found mutational profiles that were different between these patient groups. The differences concentrated among genes with low mutation frequencies in breast cancer. SIGNIFICANCE: The study describes differences in tumor mutational profiles between African American, European American, and Kenyan breast cancer patients. It also investigates how these profiles may relate to the tumor immune environment and the neighborhood environment in which the patients had residence. Finally, it describes an overrepresentation of ARID1A gene mutations in breast tumors of the Kenyan patients.

Estimated Prevalence, Tumor Spectrum, and Neurofibromatosis Type 1-Like Phenotype of CDKN2A-Related Melanoma-Astrocytoma Syndrome.

Importance: Knowledge about the prevalence and tumor types of CDKN2A-related melanoma-astrocytoma syndrome (MAS) is limited and could improve disease recognition. Objective: To estimate the prevalence and describe the tumor types of MAS. Design, Setting, and Participants: This retrospective cohort study analyzed all available MAS cases from medical centers in the US (2 sites) and Europe (2 sites) and from biomedical population genomic databases (UK Biobank [United Kingdom], Geisinger MyCode [US]) between January 1, 1976, and December 31, 2020. Patients with MAS with CDKN2A germline pathogenic variants and 1 or more neural tumors were included. Data were analyzed from June 1, 2022, to January 31, 2023. Main Outcomes and Measures: Disease prevalence and tumor frequency. Results: Prevalence of MAS ranged from 1 in 170 503 (n = 1 case; 95% CI, 1:30 098-1:965 887) in Geisinger MyCode (n = 170 503; mean [SD] age, 58.9 [19.1] years; 60.6% women; 96.2% White) to 1 in 39 149 (n = 12 cases; 95% CI, 1:22 396-1:68 434) in UK Biobank (n = 469 789; mean [SD] age, 70.0 [8.0] years; 54.2% women; 94.8% White). Among UK Biobank patients with MAS (n = 12) identified using an unbiased genomic ascertainment approach, brain neoplasms (4 of 12, 33%; 1 glioblastoma, 1 gliosarcoma, 1 astrocytoma, 1 unspecified type) and schwannomas (3 of 12, 25%) were the most common malignant and benign neural tumors, while cutaneous melanoma (2 of 12, 17%) and head and neck squamous cell carcinoma (2 of 12, 17%) were the most common nonneural malignant neoplasms. In a separate case series of 14 patients with MAS from the US and Europe, brain neoplasms (4 of 14, 29%; 2 glioblastomas, 2 unspecified type) and malignant peripheral nerve sheath tumor (2 of 14, 14%) were the most common neural cancers, while cutaneous melanoma (4 of 14, 29%) and sarcomas (2 of 14, 14%; 1 liposarcoma, 1 unspecified type) were the most common nonneural cancers. Cutaneous neurofibromas (7 of 14, 50%) and schwannomas (2 of 14, 14%) were also common. In 1 US family, a father and son with MAS had clinical diagnoses of neurofibromatosis type 1 (NF1). Genetic testing of the son detected a pathogenic CDKN2A splicing variant (c.151-1G>C) and was negative for NF1 genetic alterations. In UK Biobank, 2 in 150 (1.3%) individuals with clinical NF1 diagnoses had likely pathogenic variants in CDKN2A, including 1 individual with no detected variants in the NF1 gene. Conclusions and Relevance: This cohort study estimates the prevalence and describes the tumors of MAS. Additional studies are needed in genetically diverse populations to further define population prevalence and disease phenotypes.

SATS: A mutational signature analyzer of targeted sequenced tumors.

Tumor mutational signatures are important in clinical decision-making and are typically analyzed using whole exome or genome sequencing (WES/WGS). However, targeted sequencing is more commonly used in clinical settings, posing challenges in mutational signature analysis due to sparse mutation data and non-overlapping targeted gene panels. We introduce SATS (Signature Analyzer for Targeted Sequencing), an analytical method that identifies mutational signatures in targeted sequenced tumors by analyzing tumor mutational burdens and accounting for different gene panels. We demonstrate through simulations and pseudo-targeted sequencing data (generated by down-sampling WES/WGS data) that SATS can accurately detect common mutational signatures with distinct profiles. Using SATS, we created a pan-cancer catalog of mutational signatures specifically tailored to targeted sequencing by analyzing 100,477 targeted sequenced tumors from the AACR Project GENIE. The catalog allows SATS to estimate signature activities even within a single sample, providing new opportunities for applying mutational signatures in clinical settings.

Association of germline variants in telomere maintenance genes (POT1, TERF2IP, ACD, and TERT) with spitzoid morphology in familial melanoma: A multi-center case series.

Background: Spitzoid morphology in familial melanoma has been associated with germline variants in POT1, a telomere maintenance gene (TMG), suggesting a link between telomere biology and spitzoid differentiation. Objective: To assess if familial melanoma cases associated with germline variants in TMG (POT1, ACD, TERF2IP, and TERT) commonly exhibit spitzoid morphology. Methods: In this case series, melanomas were classified as having spitzoid morphology if at least 3 of 4 dermatopathologists reported this finding in ≥25% of tumor cells. Logistic regression was used to calculate odds ratios (OR) of spitzoid morphology compared to familial melanomas from unmatched noncarriers that were previously reviewed by a National Cancer Institute dermatopathologist. Results: Spitzoid morphology was observed in 77% (23 of 30), 75% (3 of 4), 50% (2 of 4), and 50% (1 of 2) of melanomas from individuals with germline variants in POT1, TERF2IP, ACD, and TERT, respectively. Compared to noncarriers (n = 139 melanomas), POT1 carriers (OR = 225.1, 95% confidence interval: 51.7-980.5; P < .001) and individuals with TERF2IP, ACD, and TERT variants (OR = 82.4, 95% confidence interval: 21.3-494.6; P < .001) had increased odds of spitzoid morphology. Limitations: Findings may not be generalizable to nonfamilial melanoma cases. Conclusion: Spitzoid morphology in familial melanoma could suggest germline alteration of TMG.

DNA methylation age in paired tumor and adjacent normal breast tissue in Chinese women with breast cancer.

BACKGROUND: Few studies have examined epigenetic age acceleration (AA), the difference between DNA methylation (DNAm) predicted age and chronological age, in relation to somatic genomic features in paired cancer and normal tissue, with less work done in non-European populations. In this study, we aimed to examine DNAm age and its associations with breast cancer risk factors, subtypes, somatic genomic profiles including mutation and copy number alterations and other aging markers in breast tissue of Chinese breast cancer (BC) patients from Hong Kong. METHODS: We performed genome-wide DNA methylation profiling of 196 tumor and 188 paired adjacent normal tissue collected from Chinese BC patients in Hong Kong (HKBC) using Illumina MethylationEPIC array. The DNAm age was calculated using Horvath's pan-tissue clock model. Somatic genomic features were based on data from RNA sequencing (RNASeq), whole-exome sequencing (WES), and whole-genome sequencing (WGS). Pearson's correlation (r), Kruskal-Wallis test, and regression models were used to estimate associations of DNAm AA with somatic features and breast cancer risk factors. RESULTS: DNAm age showed a stronger correlation with chronological age in normal (Pearson r = 0.78, P < 2.2e-16) than in tumor tissue (Pearson r = 0.31, P = 7.8e-06). Although overall DNAm age or AA did not vary significantly by tissue within the same individual, luminal A tumors exhibited increased DNAm AA (P = 0.004) while HER2-enriched/basal-like tumors exhibited markedly lower DNAm AA (P = < .0001) compared with paired normal tissue. Consistent with the subtype association, tumor DNAm AA was positively correlated with ESR1 (Pearson r = 0.39, P = 6.3e-06) and PGR (Pearson r = 0.36, P = 2.4e-05) gene expression. In line with this, we found that increasing DNAm AA was associated with higher body mass index (P = 0.039) and earlier age at menarche (P = 0.035), factors that are related to cumulative exposure to estrogen. In contrast, variables indicating extensive genomic instability, such as TP53 somatic mutations, high tumor mutation/copy number alteration burden, and homologous repair deficiency were associated with lower DNAm AA. CONCLUSIONS: Our findings provide additional insights into the complexity of breast tissue aging that is associated with the interaction of hormonal, genomic, and epigenetic mechanisms in an East Asian population.

Gene Expression Profiling Identifies Two Chordoma Subtypes Associated with Distinct Molecular Mechanisms and Clinical Outcomes.

PURPOSE: Chordoma is a rare bone tumor with a high recurrence rate and limited treatment options. The aim of this study was to identify molecular subtypes of chordoma that may improve clinical management. EXPERIMENTAL DESIGN: We conducted RNA sequencing in 48 tumors from patients with Chinese skull-base chordoma and identified two major molecular subtypes. We then replicated the classification using a NanoString panel in 48 patients with chordoma from North America. RESULTS: Tumors in one subtype were more likely to have somatic mutations and reduced expression in chromatin remodeling genes, such as PBRM1 and SETD2, whereas the other subtype was characterized by the upregulation of genes in epithelial-mesenchymal transition and Sonic Hedgehog pathways. IHC staining of top differentially expressed genes between the two subtypes in 312 patients with Chinese chordoma with long-term follow-up data showed that the expression of some markers such as PTCH1 was significantly associated with survival outcomes. CONCLUSIONS: Our findings may improve the understanding of subtype-specific tumorigenesis of chordoma and inform clinical prognostication and targeted options.

Gene-Level Associations in Patients With and Without Pathogenic Germline Variants in CDKN2A and Pancreatic Cancer.

PURPOSE: Pancreatic ductal adenocarcinoma (PDAC) is a component of familial melanoma due to germline pathogenic variants (GPVs) in CDKN2A. However, it is unclear what role this gene or other genes play in its etiology. MATERIALS AND METHODS: We analyzed 189 cancer predisposition genes using parametric rare-variant association (RVA) tests and nonparametric permutation tests to identify gene-level associations in PDAC for patients with (CDKN2A+) and without (CDKN2A-) GPV. Exome sequencing was performed on 84 patients with PDAC, 47 CDKN2A+ and 37 CDKN2A-. After variant filtering, various RVA tests and permutation tests were run separately by CDKN2A status. Genes with the strongest nominal associations were evaluated in patients with PDAC from The Cancer Genome Atlas and the UK Biobank (UKB). A secondary analysis including only GPV from UKB was also performed. RESULTS: In RVA tests, ERCC4 and RET showed the most compelling evidence as plausible PDAC candidate genes for CDKN2A+ patients. In contrast, the findings in CDKN2A- patients provided evidence for HMBS, EPCAM, and MRE11 as potential new candidate genes and confirmed ATM, BRCA2, and PALB2 as PDAC genes, consistent with findings in The Cancer Genome Atlas and the UKB. As expected, CDKN2A- patients were more likely to harbor GPVs from the 189 genes investigated. When including only GPVs from UKB, significant associations with PDAC were seen for ATM, BRCA2, and CDKN2A. CONCLUSION: These results suggest that variants in other genes likely play a role in PDAC in all patients and that PDAC in CDKN2A+ patients has a distinct etiology from PDAC in CDKN2A- patients.

Impact of Transcript (p16/p14ARF) Alteration on Cancer Risk in CDKN2A Germline Pathogenic Variant Carriers.

BACKGROUND: Few studies have evaluated the relationship between CDKN2A germline pathogenic variants (GPV), transcript (p16/p14ARF) alteration, and cancer risk. METHODS: Standardized incidence ratios (SIRs) comparing cancer risk with the general population were calculated for 385 CDKN2A GPV carriers from 2 large cohorts (259 United States and 126 Swedish individuals) using Poisson regression; statistical significance was defined as P less than .002 (Bonferroni correction). Cumulative incidence is reported for melanoma and nonmelanoma cancer. RESULTS: Incidence was increased for melanoma (SIR = 159.8, 95% confidence interval [CI] = 132.1 to 193.2), pancreatic cancer (SIR = 24.1, 95% CI = 14.7 to 39.4), head and neck squamous cell carcinoma (SIR = 16.2, 95% CI = 9.5 to 27.6), and lung cancer (SIR = 5.6, 95% CI = 3.4 to 9.1) in GPV carriers. Similar associations were observed with p16 alteration. Combined p16 and p14ARF alteration was associated with increased incidence of esophageal cancer (SIR = 16.7, 95% CI = 5.7 to 48.9) and malignant peripheral nerve sheath tumor (SIR = 113.0, 95% CI = 16.4 to 780.9), although cancer events were limited (n < 5 for each malignancy). Cumulative incidence at age 70 years for melanoma and nonmelanoma cancer was 68.3% (95% CI = 68.0% to 68.6%) and 35.2% (95% CI = 34.9% to 35.6%), respectively. A total 89% of smoking-related cancers (lung, head and neck squamous cell carcinoma, pancreatic, esophageal) occurred in ever smokers. CONCLUSION: These findings highlight the impact of p16 and p14ARF alteration on cancer risk. Smoking was an important risk factor for smoking-related cancers in our study.

Rare germline variants in PALB2 and BRCA2 in familial and sporadic chordoma.

Chordoma is a rare bone tumor with genetic risk factors largely unknown. We conducted a whole-exome sequencing (WES) analysis of germline DNA from 19 familial chordoma cases in five pedigrees and 137 sporadic chordoma patients and identified 17 rare germline variants in PALB2 and BRCA2, whose products play essential roles in homologous recombination (HR) and tumor suppression. One PALB2 variant showed disease cosegregation in a family with four affected people or obligate gene carrier. Chordoma cases had a significantly increased burden of rare variants in both genes when compared to population-based controls. Four of the six PALB2 variants identified from chordoma patients modestly affected HR function and three of the 11 BRCA2 variants caused loss of function in experimental assays. These results, together with previous reports of abnormal morphology and Brachyury expression of the notochord in Palb2 knockout mouse embryos and genomic signatures associated with HR defect and HR gene mutations in advanced chordomas, suggest that germline mutations in PALB2 and BRCA2 may increase chordoma susceptibility. Our data shed light on the etiology of chordoma and support the previous finding that PARP-1 inhibitors may be a potential therapy for some chordoma patients.

Rare germline deleterious variants increase susceptibility for lung cancer.

Although multiple common susceptibility loci for lung cancer (LC) have been identified by genome-wide association studies, they can explain only a small portion of heritability. The etiological contribution of rare deleterious variants (RDVs) to LC risk is not fully characterized and may account for part of the missing heritability. Here, we sequenced the whole exomes of 2777 participants from the Environment and Genetics in Lung cancer Etiology study, a homogenous population including 1461 LC cases and 1316 controls. In single-variant analyses, we identified a new RDV, rs77187983 [EHBP1, odds ratio (OR) = 3.13, 95% confidence interval (CI) = 1.34-7.30, P = 0.008] and replicated two previously reported RDVs, rs11571833 (BRCA2, OR = 2.18; 95% CI = 1.25-3.81, P = 0.006) and rs752672077 (MPZL2, OR = 3.70, 95% CI = 1.04-13.15, P = 0.044). In gene-based analyses, we confirmed BRCA2 (P = 0.007) and ATM (P = 0.014) associations with LC risk and identified TRIB3 (P = 0.009), involved in maintaining genome stability and DNA repair, as a new candidate susceptibility gene. Furthermore, cases were enriched with RDVs in homologous recombination repair [carrier frequency (CF) = 22.9% versus 19.5%, P = 0.017] and Fanconi anemia (CF = 12.5% versus 10.2%, P = 0.036) pathways. Our results were not significant after multiple testing corrections but were enriched in cases versus controls from large scale public biobank resources, including The Cancer Genome Atlas, FinnGen and UK Biobank. Our study identifies novel candidate genes and highlights the importance of RDVs in DNA repair-related genes for LC susceptibility. These findings improve our understanding of LC heritability and may contribute to the development of risk stratification and prevention strategies.

SUITOR: Selecting the number of mutational signatures through cross-validation.

For de novo mutational signature analysis, the critical first step is to decide how many signatures should be expected in a cancer genomics study. An incorrect number could mislead downstream analyses. Here we present SUITOR (Selecting the nUmber of mutatIonal signaTures thrOugh cRoss-validation), an unsupervised cross-validation method that requires little assumptions and no numerical approximations to select the optimal number of signatures without overfitting the data. In vitro studies and in silico simulations demonstrated that SUITOR can correctly identify signatures, some of which were missed by other widely used methods. Applied to 2,540 whole-genome sequenced tumors across 22 cancer types, SUITOR selected signatures with the smallest prediction errors and almost all signatures of breast cancer selected by SUITOR were validated in an independent breast cancer study. SUITOR is a powerful tool to select the optimal number of mutational signatures, facilitating downstream analyses with etiological or therapeutic importance.

Association of Genetic Ancestry With Terminal Duct Lobular Unit Involution Among Healthy Women.

Reduced age-related terminal duct lobular unit (TDLU) involution has been linked to increased breast cancer risk and triple-negative breast cancer. Associations of TDLU involution levels with race and ethnicity remain incompletely explored. Herein, we examined the association between genetic ancestry and TDLU involution in normal breast tissue donated by 2014 healthy women in the United States. Women of African ancestry were more likely than European women to have increased TDLU counts (odds ratio [OR]trend = 1.36, 95% confidence interval [CI] = 1.07 to 1.74), acini counts per TDLU (OR = 1.47, 95% CI = 1.06 to 2.03), and median TDLU span (ORtrend = 1.44, 95% CI = 1.08 to 1.91), indicating lower involution, whereas East Asian descendants were associated with decreased TDLU counts (ORtrend = 0.52, 95% CI = 0.35 to 0.78) after controlling for potential confounders. These associations are consistent with the racial variations in incidence rates of triple-negative breast cancer in the United States and suggest opportunities for future work examining whether TDLU involution may mediate the racial differences in subtype-specific breast cancer risk.

Integrated Analysis of Coexpression and Exome Sequencing to Prioritize Susceptibility Genes for Familial Cutaneous Melanoma.

The application of whole-exome sequencing has led to the identification of high- and moderate-risk variants that contribute to cutaneous melanoma susceptibility. However, confirming disease-causing variants remains challenging. We applied a gene coexpression network analysis to prioritize the candidate genes identified from whole-exome sequencing of 34 melanoma-prone families, with at least three affected members sequenced per family (N = 119 cases). A coexpression network was constructed from genotype-tissue expression project, skin melanoma from the cancer genome atlas, and primary melanocyte cultures. We performed module-specific enrichment and focused on modules associated with pigmentation processes because they are the best-studied and most well-known risk factors for melanoma susceptibility. We found that pigmentation-associated modules across the four expression datasets examined were enriched for well-known melanoma susceptibility genes plus genes associated with pigmentation. We also used network properties to prioritize genes within pigmentation modules as candidate susceptibility genes. Integrating information from coexpression network analysis and variant prioritization, we identified 36 genes (such as DCT, TPCN2, TRPM1, ATP10A, and EPHA5) as potential melanoma risk genes in the families. Our approach also allowed us to link families with private gene mutations on the basis of gene coexpression patterns and thereby may provide an innovative perspective in gene identification in high-risk families.

Comparison of somatic mutation landscapes in Chinese versus European breast cancer patients.

Recent genomic studies suggest that Asian breast cancer (BC) may have distinct somatic features; however, most comparisons of BC genomic features across populations did not account for differences in age, subtype, and sequencing methods. In this study, we analyzed whole-exome sequencing (WES) data to characterize somatic copy number alterations (SCNAs) and mutation profiles in 98 Hong Kong BC (HKBC) patients and compared with those from The Cancer Genome Atlas of European ancestry (TCGA-EA, N = 686), which had similar distributions of age at diagnosis and PAM50 subtypes as in HKBC. We developed a two-sample Poisson model to compare driver gene selection pressure, which reflects the effect sizes of cancer driver genes, while accounting for differences in sample size, sequencing platforms, depths, and mutation calling methods. We found that somatic mutation and SCNA profiles were overall very similar between HKBC and TCGA-EA. The selection pressure for small insertions and deletions (indels) in GATA3 (false discovery rate (FDR) corrected p < 0.01) and single-nucleotide variants (SNVs) in TP53 (nominal p = 0.02, FDR corrected p = 0.28) was lower in HKBC than in TCGA-EA. Among the 13 signatures of single-base substitutions (SBS) that are common in BC, we found a suggestively higher contribution of SBS18 and a lower contribution of SBS1 in HKBC than in TCGA-EA, while the two APOBEC-induced signatures showed similar prevalence. Our results suggest that the genomic landscape of BC was largely very similar between HKBC and TCGA-EA, despite suggestive differences in some driver genes and mutational signatures that warrant future investigations in large and diverse Asian populations.

Novel loss-of-function variant in DENND5A impedes melanosomal cargo transport and predisposes to familial cutaneous melanoma.

PURPOSE: More than half of the familial cutaneous melanomas have unknown genetic predisposition. This study aims at characterizing a novel melanoma susceptibility gene. METHODS: We performed exome and targeted sequencing in melanoma-prone families without any known melanoma susceptibility genes. We analyzed the expression of candidate gene DENND5A in melanoma samples in relation to pigmentation and UV signature. Functional studies were carried out using microscopic approaches and zebrafish model. RESULTS: We identified a novel DENND5A truncating variant that segregated with melanoma in a Swedish family and 2 additional rare DENND5A variants, 1 of which segregated with the disease in an American family. We found that DENND5A is significantly enriched in pigmented melanoma tissue. Our functional studies show that loss of DENND5A function leads to decrease in melanin content in vitro and pigmentation defects in vivo. Mechanistically, harboring the truncating variant or being suppressed leads to DENND5A losing its interaction with SNX1 and its ability to transport the SNX1-associated vesicles from melanosomes. Consequently, untethered SNX1-premelanosome protein and redundant tyrosinase are redirected to lysosomal degradation by default, causing decrease in melanin content. CONCLUSION: Our findings provide evidence of a physiological role of DENND5A in the skin context and link its variants to melanoma susceptibility.

Prognostic significance of tumor-infiltrating lymphocytes in premenopausal, luminal breast cancer treated with adjuvant endocrine therapy.

PURPOSE: Tumor-infiltrating lymphocytes (TILs) have strong prognostic value in triple-negative and HER2-enriched breast cancer, but their prognostic role in luminal breast cancer (LBC) is less clear. Here, we assessed the overall TIL levels and CD8+ T-cells in relation to the prognosis of LBC patients from China. METHODS AND RESULTS: A total of 596 patients with LBC who were premenopausal and treated with adjuvant endocrine therapy were included. Among them, 160 cases were evaluated for CD8 by immunohistochemical (IHC) staining. Whole-section hematoxylin and eosin and IHC staining were visually assessed for stromal TILs (sTILs), stromal CD8+ T-cells (sCD8), and intratumoral CD8+ T-cells (iCD8). Multivariable analyses were used to test the associations between TILs and disease-free survival (DFS) and overall survival (OS) with the adjustment for clinicopathologic characteristics and treatment. High sTILs (≥10%) were associated with high histologic grade (P<0.001), luminal B/HER2- (P<0.001), luminal B/HER2+ subtype (P=0.002), and high Ki67 expression (≥25%; P=0.014). Similar associations were observed for sCD8 but not for iCD8. While sTILs and sCD8 were not associated with either DFS or OS, the presence of iCD8 (≥1%) was associated with better DFS in both univariate (HR=0.51, 95% CI 0.26-0.96, P=0.042) and multivariate (HR=0.48, 95% CI 0.25-0.92, P=0.027) analyses. Similar but less significant associations were found for iCD8 and OS (adjusted HR=0.35, 95% CI 0.11-1.10, P=0.073). CONCLUSIONS: Among Chinese premenopausal patients with LBC, iCD8 demonstrated suggestive associations with favorable outcome. In contrast, although sTILs and sCD8 were associated with more aggressive tumor features, they did not appear to be associated with clinical outcome.

Quantitative Mammographic Density Measurements and Molecular Subtypes in Chinese Women With Breast Cancer.

Background: Studies investigating associations between mammographic density (MD) and breast cancer subtypes have generated mixed results. We previously showed that having extremely dense breasts was associated with the human epidermal growth factor receptor-2 (HER2)-enriched subtype in Chinese breast cancer patients. Methods: In this study, we reevaluated the MD-subtype association in 1549 Chinese breast cancer patients, using VolparaDensity software to obtain quantitative MD measures. All statistical tests were 2-sided. Results: Compared with women with luminal A tumors, women with luminal B/HER2- (odds ratio [OR] = 1.20, 95% confidence interval [CI] = 1.04 to 1.38; P = .01), luminal B/HER2+ (OR = 1.22, 95% CI = 1.03 to 1.46; P = .03), and HER2-enriched tumors (OR = 1.30, 95% CI = 1.06 to 1.59; P = .01) had higher fibroglandular dense volume. These associations were stronger in patients with smaller tumors (<2 cm). In contrast, the triple-negative subtype was associated with lower nondense volume (OR = 0.82, 95% CI = 0.68 to 0.99; P = .04), and the association was only seen among older women (age 50 years or older). Conclusion: Although biological mechanisms remain to be investigated, the associations for the HER2-enriched and luminal B subtypes with increasing MD may partially explain the higher prevalence of luminal B and HER2+ breast cancers previously reported in Asian women.

Impact of breast cancer risk factors on clinically relevant prognostic biomarkers for primary breast cancer.

PURPOSE: In addition to impacting incidence, risk factors for breast cancer may also influence recurrence and survival from the disease. However, it is unclear how these factors affect combinatorial biomarkers for aiding treatment decision-making in breast cancer. METHODS: Patients were 8179 women with histologically confirmed invasive breast cancer, diagnosed and treated in a large cancer hospital in Beijing, China. Individual clinicopathological (tumor size, grade, lymph nodes) and immunohistochemical (IHC: ER, PR, HER2, KI67) markers were used to define clinically relevant combinatorial prognostic biomarkers, including the Nottingham Prognostic Index (NPI: combining size, grade, nodes) and IHC4 score (combining ER, PR, HER2, KI67). Odds ratios (ORs) and 95% confidence intervals (CIs) for associations between breast cancer risk factors and quartiles (Q1-Q4) of NPI and IHC4 were assessed in multivariable polytomous logistic regression models. RESULTS: Overall, increasing parity (ORtrend(95% CI) = 1.20(1.05-1.37);Ptrend = 0.007), overweight (OR(95% CI)vs normal = 1.60(1.29-1.98)), and obesity (OR(95% CI) vs normal = 2.12(1.43-3.14)) were associated with higher likelihood of developing tumors with high (Q4) versus low (Q1) NPI score. Conversely, increasing age (ORtrend(95% CI) = 0.75(0.66-0.84);Ptrend < 0.001) and positive family history of breast cancer (FHBC) (OR(95% CI) = 0.66(0.45-0.95)) were inversely associated with NPI. Only body mass index (BMI) was associated with IHC4, with overweight (OR(95% CI) vs normal = 0.82(0.66-1.02)) and obese (OR(95% CI) vs normal = 0.52(0.36-0.76)) women less likely to develop high IHC4 tumors. Notably, elevated BMI was associated with higher NPI irrespective of hormone receptor-expression status. CONCLUSIONS: Our findings indicate that factors affecting breast cancer incidence, particularly age, parity, FHBC, and BMI, may impact clinically relevant prognostic biomarkers with implications for surveillance, prognostication, and counseling.

Breast cancer risk factors in relation to molecular subtypes in breast cancer patients from Kenya.

BACKGROUND: Few studies have investigated risk factor heterogeneity by molecular subtypes in indigenous African populations where prevalence of traditional breast cancer (BC) risk factors, genetic background, and environmental exposures show marked differences compared to European ancestry populations. METHODS: We conducted a case-only analysis of 838 pathologically confirmed BC cases recruited from 5 groups of public, faith-based, and private institutions across Kenya between March 2012 to May 2015. Centralized pathology review and immunohistochemistry (IHC) for key markers (ER, PR, HER2, EGFR, CK5-6, and Ki67) was performed to define subtypes. Risk factor data was collected at time of diagnosis through a questionnaire. Multivariable polytomous logistic regression models were used to determine associations between BC risk factors and tumor molecular subtypes, adjusted for clinical characteristics and risk factors. RESULTS: The median age at menarche and first pregnancy were 14 and 21 years, median number of children was 3, and breastfeeding duration was 62 months per child. Distribution of molecular subtypes for luminal A, luminal B, HER2-enriched, and triple negative (TN) breast cancers was 34.8%, 35.8%, 10.7%, and 18.6%, respectively. After adjusting for covariates, compared to patients with ER-positive tumors, ER-negative patients were more likely to have higher parity (OR = 2.03, 95% CI = (1.11, 3.72), p = 0.021, comparing ≥ 5 to ≤ 2 children). Compared to patients with luminal A tumors, luminal B patients were more likely to have lower parity (OR = 0.45, 95% CI = 0.23, 0.87, p = 0.018, comparing ≥ 5 to ≤ 2 children); HER2-enriched patients were less likely to be obese (OR = 0.36, 95% CI = 0.16, 0.81, p = 0.013) or older age at menopause (OR = 0.38, 95% CI = 0.15, 0.997, p = 0.049). Body mass index (BMI), either overall or by menopausal status, did not vary significantly by ER status. Overall, cumulative or average breastfeeding duration did not vary significantly across subtypes. CONCLUSIONS: In Kenya, we found associations between parity-related risk factors and ER status consistent with observations in European ancestry populations, but differing associations with BMI and breastfeeding. Inclusion of diverse populations in cancer etiology studies is needed to develop population and subtype-specific risk prediction/prevention strategies.

Rare Germline Variants in Chordoma-Related Genes and Chordoma Susceptibility.

BACKGROUND: Chordoma is a rare bone cancer with an unknown etiology. TBXT is the only chordoma susceptibility gene identified to date; germline single nucleotide variants and copy number variants in TBXT have been associated with chordoma susceptibility in familial and sporadic chordoma. However, the genetic susceptibility of chordoma remains largely unknown. In this study, we investigated rare germline genetic variants in genes involved in TBXT/chordoma-related signaling pathways and other biological processes in chordoma patients from North America and China. METHODS: We identified variants that were very rare in general population and internal control datasets and showed evidence for pathogenicity in 265 genes in a whole exome sequencing (WES) dataset of 138 chordoma patients of European ancestry and in a whole genome sequencing (WGS) dataset of 80 Chinese patients with skull base chordoma. RESULTS: Rare and likely pathogenic variants were identified in 32 of 138 European ancestry patients (23%), including genes that are part of notochord development, PI3K/AKT/mTOR, Sonic Hedgehog, SWI/SNF complex and mesoderm development pathways. Rare pathogenic variants in COL2A1, EXT1, PDK1, LRP2, TBXT and TSC2, among others, were also observed in Chinese patients. CONCLUSION: We identified several rare loss-of-function and predicted deleterious missense variants in germline DNA from patients with chordoma, which may influence chordoma predisposition and reflect a complex susceptibility, warranting further investigation in large studies.